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SUMMARY 


The photon spectrum at 29-30 km above sea-level and 64°-65° geomagnetic latitude has 
been determined by means of the nuclear emulsion technique, the multiple scattering of the 
electron-pairs, produced when the photons materialize in the emulsion being measured. Distortion, 
noise, energy loss of the electrons, outscattering and materialization probability have been taken 
into account. The relation between the registered photon spectrum and the z°-meson spectrum 
has been studied. The exponent obtained for the differential energy spectrum of the 7°-mesons is 
1.5 + 0.2 in the total energy interval from 0.25 GeV to 3 GeV and 2.8 + 0.4 for energies greater 
than 3 GeV. The last exponent may suffer from systematic errors due to distortion, spurious scat- 
tering and scanning efficiency, which are difficult to estimate. 


A. Some theoretical relations between the 2°-meson and the 
photon spectra 


I. Introduction 


If photographic emulsions are exposed at high altitude to cosmic rays, there will 
appear in them after development a great number of electron-pairs, which are almost 
certainly created by the photons resulting from the decay of the z°-mesons, either 
directly or through an intervening electron-photon cascade. The relation between 
the z°-meson production energy spectrum and the photon spectrum resulting from 
the decay of the z°-mesons has been studied by Carlson e¢ al. (1), Oehme (2) and 
Tiapkin (3) among others. In this part of the present work we will deduce a rigorous 
mathematical relation which makes it possible to derive the energy spectrum of the 
z®-mesons from the distribution in zenith angle and energy of the pairs observed in 
a certain volume of the photographic emulsion. 


II. The general relation 


We start with the deduction of the photon spectrum N(£, 6,) at the point Q in 
fig. 1 from the z°-meson production spectra S(¢, 0,7, 0,) at the points P = (r, 6, Py): 
N(E,0,)dEH dw, is the number of photons per second and em” with the energy 
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Fig. 1. Notations of the angles used when considering the relation between the photon-spectrum 
at Q and the z°-meson production spectrum at P. 


in (Z, E +d £) and the zenith angle within the solid angle dw, around 6, and S(e, 0, 
r,0,)dV de dq is the number of °-mesons produced per second in the volume dV 
at the point P with the total energy in (¢,¢+de) and the zenith angle within dw 
around 9. Let H (L’, r, 6,) dV d E’dw be the number of photons produced by decaying 
z°-mesons per second in the volume dV at the point P with the energy in (H’, H’ + 
dE’) and the direction towards Q within the solid angle dw. The relation between 1 
and S has been published without derivation by Tiapkin (3). Because its validity does 
not seem to be self-evident, a derivation is outlined here, the notations introduced 
in Fig. 1 being used. A photon with the energy H’ can be created by a meson with 
the energy ¢ if the angle @ between their directions satisfies the relation 


e— 4/2 E’ ; 
Vea (1) 


where é is the rest energy of the °-meson. This relation is easily obtained, if one 
considers that the two decay photons have the same energy } ¢,) and opposite direc- _ 
tions in the center of mass system. The probability {(H’|e)dH’ that a meson with 


the energy ¢ decays with one photon having its energy in the interval (H’, E’ + dE’) 
has the simple form 


cos } = 


{(E'|e)= 2/Ver—e5 if 3 (e—Ve®— 0) < B' < 3 (e + Ve*— a8) 
0 in other cases. 


(2) 


The probability that one photon will appear in the interval (4 (e — Ve £), $€) is 
one. A photon in this energy interval is always accompanied by a photon in the 
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interval (4¢, }(¢ +V 2 — €0)). From this probability distribution follows that a photon 
with the energy H’ can result from the decay of a meson with an energy greater than 


2 
€0 


4 Kf 


e(H’) =H’ + (3) 


and only from such a meson. The relations (2) and (3) are derived in ref. (1). Let 
f(E’, 6,|¢, 0, ~) dE’d0, be the probability that one of the photons, emitted in the 
decay of a x°-meson with the total energy ¢ and the direction (6, gy), will have its 
energy 4’ within (H’, H’ + dE’) and its zenith angle within (0,, 0, +d6,). Because 
the distance travelled by a z°-meson before decay is negligible, the relation between 
H(E’,r,6,) and S(e,6,7,6,) is 


H(E',r,0,)dV dE’ 2xsin 0, d0, 
= { S(e,0,r,0,) dV deda f(E", 0, |e, 0, ~) dE’ dd, 


€,0,2 


where the integration should be performed over those ¢, 9 and y which can contribute. 
The expression can be simplified to 


i IY 
H(E’,r, vs) oor tit | S(e, 0,7, 0,) f (E’, 0,| e, 0, p)dedo. (4) 
€,0,P 


We can write 
f(H’, 0, |¢, 0, p) =f(E’ |e) (0,| EZ’, ¢, 8) (5) 


where f(H’|e) has been defined above (2) and where /(6,| H’, ¢,0) d0, is the probabil- 
ity that a photon with the energy HZ’, which is emitted from a decaying z°-meson 
with the energy ¢ and the zenith angle 0, has its zenith angle within (6,, 0, +d6,). 
For EH’, ¢ and @ given 6, is determined from the relation 


cos 0, = cos} cos§+sin# sinOcosB (6) 


where # is a function of H’ and « given by the relation (1). The probability g(f| Z’, 
€,0) df that f will take a value within (6, 6 +d) for fixed values of H’, ¢ and @ is 
d B/x. Because under these conditions 6, is a function of 6 only, we have 


{(0,,|H’, ¢, 0) d0, =g9(B|£’, ¢, 0) dB 
which yields 
{(0,|E’, e, 0) = 


ale 


Ce 
do, 
or after differentiation of (6) and elimination of df/d6,: 


sin 6, 
az sin} sin # sin B 


} (0,|H’, €, 0) = 


From (2), (3), (4), (5) and (7) we obtain 
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1 f de [S(e,0,7,0)) do 
H (E’,r, 9,) =a V2—e2 J sind sin 6 sin p 
0 


e(E’) & 


Because S is assumed to be independent of ¢, this angle can be integrated out, and 
putting dw =sin0d0 dq we obtain 


20 de Sd 
, he . 8 
Hitt, By) cA | el aes eee is 


e(E’) 6 


There now remains a change of variables. We want to replace 6 by y and integrate 
over the latter variable instead of over 0. For this purpose we can use the spherical 
trigonometric relation 

cos 0 = cos # cos 6, + sin & sin 6, cos p (8a) 


which after differentiation gives 
sin 9 d6 =sin # sin 0, sin yp dy. 
If we make use of another relation from spherical trigonometry, 
sin 6 sin 6 = sin 6, sin p 
we obtain 


fp Lake 
?™ sind sin B 


When this is substituted into (8), we obtain the integration limit upon geometrical 
consideration of fig. 1 and 


4 


: 2 pod 
H(B',r,8,) == | a] Seer aay (9) 
™~ €0 
€(E’) 0 


which is the final expression for the relation between H and S. In order to get N (E, 0,), 
let us consider the photons that pass through the area dS in fig. 2 with the energy 
in (HZ, # + d£) and the direction within the solid angle dw, around (6,, y,) and orig- 


inate from decay photons in the distance-interval (7,7 +dr). If the lateral spread 


of the shower is left out of consideration, the later photons must have been created 
in the volume 


dV =rdw,dr 
and the number per second of the former photons is 
dN (E,0,)dHdw,dS=dE { dB’ y(E, E',r,0,) H(E’,r,0,)dV dow 
E’=E 
where dw =dS/r? and where y(E, E’,r, 6,) dE is the average number of photons 


with the energy in the interval (H, H + dE) that arrive at Q as a result of a photon 
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ZENITH 


Fig. 2. The solid angles occurring in the deduction of the 
relation between the photon-spectrum at Q and the 7°-meson 
production spectrum in the vicinity. 


with the energy H’ starting at P moving towards Q. y can be obtained from shower 
theory. The above expression may be simplified to 


dN (E,0,)=dr | dE’ y (HE, E’,r, 0,) H (E’, r, 0,) 
E 
and from this we obtain 


N (E, 6,)= [dr {dE y (E, E’,r, 0,) H(E’, 7, 0,) (10) 
0 E 


which is the final expression for the relation between N and H. 

The photons are detected by inspection of the pairs that they produce in a photo- 
graphic plate. The finite extension of the detector may be taken into account in the 
following way. Let N(Q, H, 6,) be the photon spectrum at the point Q of the detector 
and let P(Q, LH, 0,) dV dE dw, be the number of pairs that are created per unit time 
in the volume dV around Q with the energy in (H,H#+d£) and the zenith angle 
within the solid angle dw, around 6,. The direction of a pair is supposed to be the 
same as that of the parent photon. Let o(#) be the cross-section per atom for pair- 
production and » the number of atoms per cm’. We have 


P(Q, EL, 6,) =n oN (Q, E, 6,). (11) 


The number of pairs per unit energy and solid angle created in the volume V during 
the time 7’ is 
R(B,0,)= { P(Q, H,6,) dV dt (12) 
V,T 


In principle this relation between the pair spectrum in a volume and the z°-meson 
production spectrum can be used when deducing the latter from the former. For 


instance one can assume 
S =a,8, a5 AyS5 Se ole 


where S,, S,,... are suitably chosen meson production spectra and aj, ag, ... are 
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constants. Because of the linearity of the relations (12), (10) and (9) the pair spectrum 
will take the form 

where R, is the pair spectrum corresponding to S;. The coefficients a; can be deter- 
mined from the experimental pair spectrum. 


III. Some special cases 


Finally we want to consider the relation between the photon spectrum N (£, 6,) 
at one point and the z°-meson production spectrum S(e,0,7,6,) in the vicinity. 
We deal with some special cases, neglecting cascade processes. 

If the z°-meson production spectrum is isotropic and has an energy dependence, 
which is independent of the position of the point P = (r, 6,, ~,), we can put 


(6; 0:9',.05)= 1 (F053) FOL ey 


Because S is now independent of y for fixed values of ¢, H and 0, (10) and (9) take 
the form 


” P(e) de 
N (EL, 6,) =k (6,) a 13 
3 : Ve*— & 
€(£) 
where k(6,) =2 f f(r, 0,) dr. 
i) 
From this the well known relationship 
2 
€0 ] aN 
F\E = 
( +7) wey * Sz a4 


can be deduced [ref. (1)]. 


; The following simple model seems to be useful in the study of the cosmic radiation 
in the upper part of the atmosphere: 


{f(r 6,)F(e) if 0<4x 


Cae 
lo if O>42. 
If we put I (e, r, 6,) = iKic 6,7, 0,) dy, 
we obtain: : 
1 if 0<0,<3n2, 0<8<i4n-6, 
Yo (9, 9), * » $a—-O,<0<4240, 
I=F-] 0 a . » $2+0,<0<an 
0 » 30<0,<2, 0<8<0,-42 
Yo(I,9) 5, 4, sO, $m <<32/2-6,. 
a es 5 , 32/2-0,<0<a2 


where wo(#, 0,) = arccos(— cot 6, cot ?). 
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In order to obtain these results we have to find out, for what values of wy we have 
0 <6 < $a. This depends on 6, and #, where # is a function of H’ and ¢ given by formula 
(1). From fig. 1 we see that 9 = Onin = |, —#0| when y =0. When y increases to z, 
§ increases to Omax = 0, +9 or 22 — (6, +@), the expression, which is:smaller than z, 
being correct. From (8a) the condition for 6 = 42 is obtained, namely 


cos # cos 6, + sin sin 6, cos yp = 0, 
which can be written as 
Y = Wo(V, 8,) = arccos ( — cot # cot 6,). 


When |cot # cot 6,| >1, y does not exist and all values of y between 0 and z 
can contribute. If on the other hand Oni, > 4.2, there are no values of y that can con- 
tribute. The result written down above is obtained as an immediate consequence of 
these two conditions. It should be noted that (13) and (14) hold in the two cases, 
when [=F fz. 


B. An experimental investigation of the photon and <°-meson 
spectra at high altitude 


I. Introduction 


The main properties of the cosmic radiation are fairly well known nowadays. 
It is most likely that the primary cosmic radiation consists of protons and heavier 
nuclei. The nucleons have a ‘differential energy spectrum N(H)d£. The exponent 


_ _ émN 
sony AE 


of this distribution is a slowly increasing function of the energy # and seems to tend 
to a constant value for great energies #. Experiments with extensive air showers 
(4) have indicated that the variation of p is very small above 14 GeV, the average 
value between 14 and 1.4- 108 GeV being 2.64 + 0.04. If, for a rough estimate, we 
assume the same exponent 2.6 between 1 GeV and infinity, we find that in this 
energy region, where multiple meson production is important, less than 10% of 
the total energy is carried by primary nucleons with an energy greater than 50 GeV. 
Only at energies greater than this the total production cross-section for nucleon- 
antinucleon-pairs, hyperons and K-mesons is believed to be of the same order of 
magnitude as the production cross-section for a-mesons (5). The available experi- 
mental determinations of the ratio of the cross-sections involve many assumptions. 
In an analysis of jets, caused by nuclear collisions of energy between 10* and 104 GeV, 
Brisbout et al. have found that only about 25% of the shower particles are not 
a-mesons (6). It can be concluded that heavy mesons, hyperons and antinucleons 
play no major part in the development of the local cosmic radiation. It is also well 
established that high energy electrons or photons give no large contribution to the 
primary cosmic radiation. As has been pointed out by Rossi (7) an experiment by 
Hulsizer shows that not more than about one percent of the primary particles with 
energies above 4.5 GeV are electrons or photons. The electron-photon component 
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in the local cosmic radiation is instead considered to come mainly from the decay of 
7°-mesons, which in high energy nuclear reactions, where multiple meson production 
is of importance, have been found to contribute * /, of the total number of produced 
z-mesons (1, 8). There are other established reactions, which can give rise to high 
energy photons or electrons, for instance: 


192° Ao sey 

2) we Btyty 

3) Inner Bremsstrahlung (Abstrahlung) 
4) Bremsstrahlung 


The first three of these reactions are certainly not very important at high altitude. 
The cross-section for the creation of 2° is certainly much smaller than the cross-section 
for the production of z°-mesons in the most important energy region. The reaction 
2) is certainly not important at high altitude, among other things because of the 
comparatively long lifetime of the 4-meson. The energy contribution to the photon- 
component from inner Bremsstrahlung, emitted when the charged z-mesons are 
created is less than 2 % of the contribution from the decay of z°-mesons (2, 9). Brems- 
strahlung from particles heavier than electrons is negligible. Bremsstrahlung from 
electrons cannot be ruled out at once, being an important part of the electron-photon 
cascade process, the possible contribution of which will be discussed in section II. 


II. Experimental methods and corrections 
1. General 


In the present investigation Ilford G5 plates, 7.5 em < 7.5 em exposed on the 15th 
of October 1950 in the United States at a geomagnetic latitude of 64°-65° have been 
used. The time spent at an altitude of more than 29 km above sea-level was 4.0 
hours and most of this time was spent between 29 km and 30 km. 

Certain areas of the plates have been scanned for electron-pairs. Two samples of 
pairs have been selected according to the projected length L of the shortest track. 
Sample I contains 56 pairs with ZL > 1 mm, found in one of the scanned areas and 
sample IT consists of all found pairs with L > 3 mm. Pairs, which go near and nearly 
parallel to a minimum track, have been rejected as they can have been created by 
Bremsstrahlung from the particle, which gave rise to the minimum track. 

The energy of the electrons has been determined from the multiple scattering. 
Fowler’s coordinate method (10) has been used when determining the mean scattering 
sagitta D*, and large scattering sagittas have been eliminated without replacement 
until none was greater than 4 D*. The measurements have been made with a type 


M 40000 microscope, manufactured by Cooke, Throughton & Simms and equipped 
with an eyepiece filar micrometer. 


2. Distortion 


The plates were not distortion-free, but because most tracks had only a small 
angle of dip, distortion was in most cases negligible. Those tracks, which have 
been measured with a cell length t > 400 jt have been corrected for C-shaped distortion 


in the following manner, though this correction has very small influence on the 
energy spectra. 
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Let D* be the measured sagitta. D* has an approximately Gaussian distribution 
with the mean value 0 and a standard deviation which we label o. Assuming a Gauss- 
ian distribution we can write: 


<Dty =0 
(Di?) =6? 


(jDth=p*-)2 « 


, | ge 
If every Dy is corrected so that = S D* =0 after the correction, the correction will 
i=1 
be too large because this quantity has a Gaussian distribution (0, o/ Vn). The 
variable 


aoe | es 


n SP 


i=1 


is distributed according to Student’s distribution with n —1 degrees of freedom, 
which distribution is asymptotically normal (0, 1). Because the distributions of 


VD? and =| D#| are similar in shape and because 


Zien Ses ieee 
<2 [DED = 2 «1D: p=n 20-22 nat /22 V<> Dy 


it is to be expected that 


> Di 
i=1 
DF 


De 


has a distribution which is asymptotically normal (0, 1) approximately. If we have NV 
tracks and not too small values of n, we thus expect 


i=1 


34 tracks have been measured with ¢t > 400 wu and the value of S?(w) was 1.59. For 
the 24 tracks which had u? <2.6 the value of S?(w) was 0.88. For the resting 10 
tracks the value of = D; was lowered according to the method described by Fowler 
(10) until w was about 1.3. After this correction the value of S?(w) was 0.98 for all 
34 tracks. 


3. Noise 


To begin with the noise was eliminated by means of measurements with two 
different cell lengths. In the elimination account was taken of the variation with 
the cell length of both the noise and the scattering constant. 
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Notations: 


p = momentum of the electron in MeV/c 
B = velocity factor of the electron 


n 
D*= >. | D;*| = measured mean sagitta 
Ni=1 

D*=<p">= <|Di | = expectation value of the measured sagitta 
D = expectation value of the true scattering sagitta 
B = expectation value of the noise scattering sagitta 
t =cell length in units of 100 uw 

1.8 D 


alan aa expectation value of the true scattering angle (D is measured in s). 
a 


pB can be determined from the well known relation 


tl 
o4 


pp=k (1) 


where k is the scattering constant. B and k are assumed to depend on ¢ in the following 
ways: 


B=B,! k=kgt? 


From (1) we obtain D= D,t??*%. 
If we assume a Gaussian distribution for D; and B;, we obtain the relation 


D2 = D*2 — B2 (2) 


For two measurements with cell lengths ¢, and ¢, resp. we then find 


Dy 7 t 2p - (3) 
prea grea (4 
ty 
An estimator D, for D, is obtained from (3) if D* is replaced by D*. Then 
B,=V Dt — D3 824 (4) 


is useful as an estimator for the noise in the small cell t,. If the noise is known to be 
B’, one can use the convenient estimator 


D2= D*2 — BB” (5) 


for the determination of D. The statistical errors and the optimum conditions have 
been discussed recently by Ekspong (11). 


The values of the scattering constant & given by Voyvodic in reference 12 have 
been used. The total energy E of the electron is related to pp in the following way: 


E=pB(1+d), a} (i +(=3) - 1) 
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CELL 
25 50 100 200 400 800 1600 yp 


Fig. 3. The result of the two noise-determinations mentioned in the text. The errors indicated are 
statistical standard errors. 


where m is the rest energy of the electron (0.51 MeV). The term d which is approxi- 
2 . 


mately equal to (F,) when this quantity is much smaller than one, could be neg- 


lected in all cases, because it is very small at high energies; for instance d = 0.0026 
when pf = 10 MeV/c. 

Measurements on a very high energy electron gave a preliminary value B’ of the 
noise B. The cells ¢, and t, were then chosen in such a way that 1’ = D’/B’ satisfied 


the following conditions: 
2 <1 for the small cell ¢, 


2<A’' <5 for the big cell ¢, 


The values of p obtained from the preliminary noise determination was used when 
estimating D, by means of the relation (3). The noise values obtained by means of 
(4) showed, however, no larger variations than might be expected statistically. 
Therefore such values of B for which 


0.15<1<1.3 (7-5) 


were selected for a determination of the noise. The result is shown in fig. 3a. A 
weighted least square fit of a straight line to the points in the log-log-plot yields 


B, = 0.1327 + 0.0008 and p = 0.216 + 0.008. 
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The errors indicated include only statistical errors. Care has been taken to keep the 
systematic and random errors as small as possible (11). During the course of the 
investigation the construction of the microscope was changed somewhat. Fig. 3b 
shows the result of a determination of the noise after this operation. In this case one 
series of values for cells between 200 and 1600 w was obtained from measurements 
on a straight line drawn on a glass-plate! in view of the fact that such a line gives no 
grain noise. The curvature of the line and the microscope stage did not need to be 
taken into account. A least square fit yielded 


B, = 0.1342 + 0.0011 w and p = 0.244 + 0.006. 


For most of the tracks D has been estimated by means of the very convenient 
formula (5) using the noise values obtained in the determinations mentioned above. 
A signal to noise ratio 1’ =3 was aimed at and 1.9</’ <6 was accepted. 


4. Energy loss 


The relatively large length of many of the tracks makes it necessary to take into 
account energy loss by Bremsstrahlung and ionization. Let us assume that measure- 
ments have been undertaken between points at the distance x, and 2, (in mm) 
from the beginning of the track (2, > a,). The correction to be added because of the 
ionization energy loss is according to Fowler (10): 


A; = 0.59 d MeV 


L_— Xy 


where d=2,+ > 


If 40 % or more of the energy of the electron seemed to get lost by Bremsstrahlung, 
the low energy part of the track was excluded. Therefore the quantity to be added 
as a correction for Bremsstrahlung was chosen as 


Too 
2 Rf (pp) 


where R(#) is a variable radiation length defined by the relation 


A=pB 


where £ is the energy of the electron and = is the average energy loss per mm 
> 


path length due to Bremsstrahlung. The correction has been chosen so that it is 
best for H = 100 MeV but the error for other occurring energies is not serious. The 
value & (co) = 30mm has been used (13), and the variation with energy of R according 
to quantum electrodynamics (14) has been taken into account. 


* The author wishes to thank Mr. L. Lundin at th i i : 
liph.quality seating-lmiee) in at the Nobel Institute in Stockholm for supplying 
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b) SAMPLE I 
€) SAMPLE I 


L>1mm 


SJ 
N 
Sa 


oy 


Gey"! 


c) ALL PAIRS IN SAMPLE I 
WITH L>S5mm 


Fig. 4. The smaller values in the histograms show the number of pairs per GeV. The greater values 
show the photon-spectrum obtained with outscattering and materialization probability taken 


into account. The smooth curve is the photon-spectrum obtained in section III, 


5. Outscattering 


The pair spectra obtained for L > 1 mm (sample I), L > 3 mm (sample II) and 


I>5mm (selected from sample II) are shown in fig. 4. They show such a large dif- 
ference in shape that there must be some systematic error. The reason for the varia- 
tion in shape must be that the probability for outscattering is much greater at low 
energy than at high energy. The probability F'(H#,x) that an electron, found to 
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F (ser) 
1.0 
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0.2 
0.1 
E 
10 20 50 100 200 500 1000 MeV 


E 
Fig. 5. The probability F (a that the track of an electron of energy HZ found to start parallel 
xVu 


to the emulsion (0.4 mm thick) attains a projected length of x or more, calculated for x = 5 mm. 


start its path in a plane parallel to the emulsion with the energy H, will travel a 
projected distance of a or more before being scattered out, has been calculated from 
the theory of Fermi (15) under the simplifying assumptions mentioned in the appen- 


E 
dix. The result is shown in fig. 5. It is proved in the appendix that F(H, x) =F (sa 
cle 


so that the calculations need be performed only for one value of x. If pairs, for 
which the projected length is longer than L for both tracks, are selected, the intensity 


2 
is considered to be reduced by outscattering with a factor F ( eel where E is the 
energy of the pair. As pointed out in the appendix, the real reduction is certainly 
somewhat greater for small energies for which F? is appreciably smaller than one. 


6. Materialization probability 


The fraction of the photons of energy HZ in a homogeneous beam that will materialize 
in a small volume dV when passing through it is proportional to o(Z) dV, where 
o(£) is the cross-section for pair creation. The photon spectrum has been obtained 


from the corresponding pair spectrum by applying to each pair with the energy H 
the statistical weight 


a(co) 
E 2 
Ae 
Fae, oe 
o(£) 


The variation of et with energy has been obtained from quantum electrody- 
namics (14). 
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Fig. 6. The distribution in projected zenith- 
angle ¢ of all measured pairs. 5 
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Ill. Experimental results 
1. The angular distribution 


The best way to analyze the result would no doubt have been according to the 
method given in section A. The statistical material does not, however, seem to be 
large enough to justify the tremendous numerical work necessary when using that 
method. Therefore a more conventional discussion of the results has been undertaken. 

The distribution in projected zenith angle ¢ of the pairs measured is shown in fig. 6. 
This angle does not differ seriously from the true zenith angle 6,. The solid angle 
dw, corresponding to an angular interval d ¢ is independent of ¢. The angular distri- 
bution differs markedly from that obtained by Carlson e¢ al. (1). The latter was iso- 
tropic for 6 < 90° and fell off rapidly to zero for ¢ > 90°. The present distribution 
shows a strongly pronounced peak at about 90°. This peak would be still more 
marked for the distribution of the number of pairs per steradian as a function of 6, 
because 6, generally has a value somewhat more near 90° than ¢ has. The difference 
between the two experimental distributions is certainly due to the fact that the 
Bristol plates were exposed at a greater atmospheric depth (21 km above sea-level) 
than the plates used in the present investigation. In the latter case the amount of 
matter in which meson-producing collisions can occur is much greater in the direction 
near 90° than for smaller zenith angles. The small central peak is believed to be due 
to another stack placed some distance above the one used in the present investiga- 
tion during the flight. Fig. 7 shows the distribution in ¢ for all scanned pairs. It is 
similar to that in fig. 6. Because most of the pairs in fig. 7 are shorter than 1 mm 
the difference between ¢ and 6, is generally greater for this distribution than for 


that in fig. 6. 
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Fig. 7. The distribution in projected zenith-angle ¢ of all scanned pairs. 


2. The energy distribution 

A photon with the energy EH can arise in the decay of 7°-mesons of energy between 
é(H#) and co, where f 
€0 


e(H)= +7 (6) 


€) is the 2°-meson rest energy. The angle # between the direction of the 7°-meson 
and the photon depends on £ and é¢ according to 


psy 
2H 

cos } = ———— (7) 
Ve? — & 


and varies from # =0 when the z°-meson energy ¢ has the value ¢(EZ) to Omax = are- 
cos 8 when ¢ = 2 # and then back to # = 0 as eco. fc is the velocity of them°-meson. 
A plot of ®max as a function of HZ is shown in fig. 8. From this we see that a photon 
with the energy 0.26 GeV or more can form an angle of at most 15° with its parent 
z-meson. Assuming that the primary cosmic radiation has an isotropic angular 
distribution in space outside the region where the earth’s influence is important, it 
seems reasonable to assume that the 2°-mesons are isotropically distributed for zenith 
angles smaller than 90°. Then the photons with energy > 0.26 GeV and zenith angle 
< 75° have all their parent z°-mesons within this isotropic distribution and if cascade 
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Fig. 8. The maximum angle #max between 
the direction of a photon of energy H# and 
its parent 7°-meson. 1 E 
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processes are neglected the differential energy-spectrum of the 2°-mesons F'(£) 
may be obtained from that of the photons N(#) from the relation 


= 7 


4H) k| dH (8) 


F(z+ 

where & is a constant. (Section A) 
Fig. 9 shows the photon spectra obtained from pairs in the projected zenith angle 
intervals (0°, 75°), (75°, 105°) and (0°, 180°). A smooth curve has been fitted to the 
(0°, 180°)-spectrum in the following way. The differential photon spectrum was 


assumed to have the form 
N(#)=A-H? 


where £ is the photon energy in GeV. Two least square fits of straight lines to the 
experimental points in the log-log-plot gave the following result: 


0.26<E<14 A=29.74+14 p=1.53+0.06 
14<E<65 A=364474 p=2.22+0.21 


The total energy represented by all pairs of energy > 3 GeV was 37 + 16 GeV. The 
statistical error includes the error due to the expected statistical fluctuation of the 
number of pairs according to Poisson’s statistics. If the same exponent p is assumed 
between 3 GeV and oo and the value from the second least square fit is used for 
N(3 GeV), we obtain 

p=2.8+0.4 for H>3 GeV. 


This result is represented by the dotted line in fig. 9c. To the statistical standard 
error given there contribute in the last case systematic errors, depending on distortion, 
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Fig. 9. Energy-spectra for photons in the projected zenith-angle intervals indicated in the dia- 

grams. The statistical standard errors are indicated, when the number of pairs in the energy- 

interval is greater than one. The straight lines and the smooth curve in diagram c) are explained 
in the text. The smooth curve in diagrams a) and b) is the same as that in diagram ec). 


spurious scattering and scanning efficiency, which are difficult to estimate objec- 
tively. Therefore the good agreement with earlier results must not be stressed too 
much in this case. After these lines had been drawn, a smooth curve was drawn 
0 


2 


in accession to them. The curve for HZ <~—° was obtained from the fact that there 
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Fig. 10. The differential z°-meson production spec- 
trum obtained from the photon-spectrum (the smooth 
curve in fig. 7c). The dots are the results of com- 
putations according to formula (8) in section B. 


01 02 OS 2 5. 10 


should be mirror symmetry in the log-log-plot about the line H -2 if the theory of 


Carlson et al. (1) is applicable. In this investigation it has not been possible to obtain 
a value of the 7°-meson mass, because the choice of L > 3 mm for most pairs makes the 
outscattering so large for the low energy part of the spectrum that its form cannot 
be determined experimentally. 

As can be seen from fig. 9 it is possible to fit the curve obtained for the (0°, 180°) 
spectrum to the two other spectra. The spectra for (0°, 75°) and (75°, 105°) do not 
differ appreciably from each other. Cascade processes in the atmosphere are negligible 
in the former interval, but may be important in the second. Because the distance 
from the edge to the center of a plate is about one radiation length cascade processes 
in the plates are never negligible, but they are considered to give the same contribu- 
tion for all directions. Since cascade processes may thus have played a greater part in 
the (75°, 105°) interval than in the (0°, 75°) interval it seems justified to conclude that 
they are not decisive for the form of the spectrum. This may be due to the fact that 
in fairly large energy intervals the photon energy spectrum obeys a power law with 
a constant exponent p. If the exponent had been the same up to infinity it would not 
have changed during the propagation of the photons through matter according to 
shower theory. In the present case the energy intervals, where the exponent may be 
considered to be constant, are probably large enough for the exponent not to change 
seriously during the one or two radiation lengths that are possible. 

In view of the above discussion it is probably permissible to use (8) when deducing 
the z°-meson production spectrum from the photon spectrum. The result is shown in 
fig. 10. The exponents obtained are 


0.25 GeV <e <3 GeV p=1.5340.2 
é2>3 GeV p=2.8+0.4 
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This result differs from the earlier ones, which have been put together by Kim (16) 
in that the transition from the smaller exponent 1.5 to the larger 2.8 does not occur 
at 1 GeV but at a larger energy 3 GeV, or perhaps even higher if the systematic errors 
mentioned above are taken into account. The abrupt transition has perhaps no real 
significance but depends on the way in which the curve for the photon spectrum has 
been drawn. The error in the exponent p seems difficult to estimate. The statistical 
errors from the corresponding straight lines fitted to the photon spectrum have 
been used above. 


APPENDIX 


Outscattering theory 


Let us introduce a Cartesian coordinate-system (x, y, v) with the xv-plane parallel 
to the boundary surface between the emulsion and the glass and situated at the 
distance z from that surface. Consider a track beginning at (0, 0, 0) in the direction 
of the x-axis (fig. 11). The probability S(H, x, y) dy that the track has its y-coordinate 
within (y, y +dy) when the 2-coordinate is x has been calculated by Fermi (14). 
The result can be written in the form 


#H ea 

e 3 
S |—, ee ey ek 1 
Ere v) V2naVx (1) 


when the scattering constant is taken to be k = 30 (grad-MeV- 


where e)= is 
- (100 u)~""), the energy ZH of the electron is measured in MeV and a and y are meas- 
ured in mm. The substitution 


W=W(—-,y) = <b (2) 
ol Oe a Var 
‘ E 
yields 8( =v) ay-f(n) aw 
a Var 
a 
where {(W)= Be 


V2x 


LASS 


Fig. 11. Illustration to the outscattering theory. 
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Fig. 12. The probability F(— +8) 04 
w Vm 

that an electron starting at the dis- 0.3 

tance z from the glass and parallel 

to it attains a projected length of x 02 

or more in an emulsion 0.4 mm thick. 

The numbers attached to the curves 9, [2105 a ee ay 

are the electron energies H in MeV, 

or enna oN oe 

0 01 0.2 0.3 0.4 mm 


The probability, F(H, x, z) that the electron is not scattered out, until it has passed 


the x-coordinate «x is 

Ue 

E E 
F Te ORCA es 8(-u)ay 

y, 
where y, = —2, Yg=t—z (=the emulsion thickness). An approximation which 
involves an underestimation of the outscattering has been made because a track 
like that in fig. 11 is considered not to be scattered out. However, the error introduced 


in this way does not seem to be very large because the electron is rarely scattered back 
into a layer which it has already left (13). The substitution (2) yields 


r(cret)- fms 


E ; 
where W,=W Be n 1=', 2. 
x Vx 


F 2) has been calculated for £=0.4 mm and «=5 mm. The result is 
cTyer 


shown in fig. 12. 

Finally let us consider the case when all possible values of z are equally probable. 
The probability that the electron will reach the z-coordinate # without being scattered 
out is 
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ryt fede 


This integration has been carried out numerically and the result is shown in fig. 5. 

The effect of energy loss, which is such that it increases the outscattering, has not 
been taken into account. In the application of the result, account has not been taken 
of the fact that the angle of dip at the beginning of the electron tracks is generally 
not zero. 
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